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! Abstract 

It is known that a plasma in a magnetic field, conceived microscop- 
(-H ■ ically as a system of point charges, can exist in a magnetized state, and 

^ f thus remain confined, inasmuch as it is in an ordered state of motion, 

with the charged particles performing gyrational motions transverse to 
the field. Here we give an estimate of a threshold, beyond which trans- 
■ verse motion become chaotic, the electrons being unable to perform 

even one gyration, so that a breakdown should occur, with complete 
loss of confinement. The estimate is obtained by the methods of pertur- 
bation theory, taking as perturbing force acting on each electron that 
^ ■ due to the so-called microfield, i.e., the electric field produced by all 

OO , the other charges. We first obtain a general relation for the threshold, 

which involves the fiuctuations of the microfield. Then, taking for such 
fluctuations the fomula given by Iglesias, Lebowitz and MacGowan for 
the model of a one component plasma with neutralizing background, 
we obtain a definite formula for the threshold, which corresponds to 
a density limit increasing as the square of the imposed magnetic field. 
Such a theoretical density limit is found to fit pretty well the empirical 
data for collapses of fusion machines. 



^ ■ 1 Introduction 



The existence of a transition from order to chaos in Hamiltonian systems, 
as a generic phenomenon occurring when a perturbation is added to an 
integrable system, is a weh established fact. This fact involves deep mathe- 
matical features (see for example |lj ) , and was made popular in the scientific 
community through the striking pictures of Henon and Heiles [2] and the 
discovery, by Izrailev and Chirikov [3], that the ordered motions found by 
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Fermi, Pasta and Ulam in their model become chaotic above a certain 
stochasticity threshold (see Figs. 4.3 and 4.5 of the review [5], or [6]). Transi- 
tions of this type were met also in the frame of plasma physics, in connection 
with the destruction of magnetic surfaces [TJ [8] , and also with the chaoticity 
thus induced on single particle motions 

On the other hand, in plasma physics a phenomenon of great relevance 
exists that is yet unexplained, and for which we propose here an explanation 
just in terms of a transition from order to chaos. We refer to the loss of 
plasma confinement, a plasma collapse that is met when the plasma density 
is increased beyond a certain density limit (see [10], Fig. 3). 

Let us recall that confinement amounts to keeping the charged particles 
away from the walls, and is actually achieved by means of a suitable magnetic 
field, the form of which depends on the concrete machine (either just a field 
imposed from outside, or a superposition of the imposed one with that due 
to a plasma current). So, when the phenomenon of destruction of magnetic 
surfaces was understood, people thought that it might play a role in explain- 
ing the breakdown occurring at the density limit (see for example ) . The 
idea is that the gyration centers should move, when the interactions are 
neglected, along the field lines, so that some peculiar topological features 
of the latter, particularly in presence of destructions of the corresponding 
magnetic surfaces, might be relevant. However, such considerations did not 
prove sufficient to explain the quick collapses of plasmas. 

Here we propose a solution of a different character, completely unrelated 
to peculiaries of the field lines, up to the point of applying even in the 
extremely idealized case in which the field is uniform, so that the field lines 
are just straight parallel lines, covering the whole space (and the plasma 
is uniform too). The phenomenon which we take into consideration is the 
existence of a magnetic pressure, which is essential in keeping the particles 
away from the walls. On the other hand, it is well known that such a pressure 
exists inasmuch as the plasma is diamagnetic, which means, in microscocpic 
terms, that each electron is equivalent to a magnetic moment, just in virtue 
of its dynamical property of performing gyrational motions transverse to the 
field lines. This is the kind of ordered motions we are referring to. 

Now, such ordered motions induced by the external magnetic field, per- 
sist indefinitely in the unperturbed case, when one neglects the perturbation 
due to the so called microfield, i.e., the microscopic electric field acting on 
each charge and due to the Coulomb interactions with all the other ones (see 
the recent review [12] )■ On the other hand the intensity of such a pertur- 
bation clearly increases with the density, and so it seems natural to expect 
that when the perturbation is large enough, i.e., at a large enough density, a 
transition to a state of chaotic motions should occur, in which diamagnetism 
is lost, together with magnetic pressure. The proposal advanced here is that 
such a kind of transition may explain the loss of confinement, at least in its 
gross features, by providing a theoretical estimate of the density limit that 
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should be compatible, as far as order of magnitude is concerned, with the 
observed ones. 

This is what we actually find out. Working with an extremely simple 
model, we predict a stochasticity threshold, which corresponds to a density 
limit that fits pretty well those observed in collapses of several kinds of fusion 
machines. 

A key point for obtaining the stochasticity threshold is the possibility of 
applying the methods of perturbation theory to macroscopic systems, as a 
plasma is. Indeed, it is well known that the classical estimates, which were 
brought to their extreme range through the work of Nekhoroshev \13\ [T^ [15] , 
don't apply in the thermodynamic limit. In the present days, however, we 
can rely on the recent progress which was obtained when the classical scheme 
was extended as to hold in a probabilistic frame |T6] , contenting oneself with 
controlling the vast majority of initial data (with respect to a given invariant 
measure), rather than all of them. Along such lines of research, perturbation 
techniques were then extended as to be applicable in the thermodynamic 
limit [I71[I81[I9]. 

Using such techniques, we show that the zeroth order of perturbation 
theory is already sufficient to produce a stochasticity threshold, which turns 
out to fit pretty well the empirical data. A characteristic feature of such 
new approaches, which combine perturbation and statistical methods, is 
that, in estimating the relaxation times or the stochasticity thresholds, the 
perturbing force enters through its fluctuations. In the present case, what 
matters is the fluctuation of the microfield, and we are in the fortunate 
situation that a theoretical formula for it is available (for the model of a one 
component plasma with neutralizing background), which was given already 
in the old work of Iglesias, Lebowitz and MacGowan [20j. This finally allows 
us to get the theoretical formula for the density limit which is proposed here. 

2 Confinement and perturbation theory 

The idea that the loss of confinement in magnetized plasmas corresponds to 
a transition from order to chaos is easily understood. Indeed an essential 
point in guaranteeing confinement is the existence of a magnetic pressure. 
Now, in a macroscopic magnetohydrodynamic description of the problem, 
the existence of a magnetic pressure is derived from the constitutive equa- 
tions of a plasma, i.e., from the assumption that the plasma be diamagnetic. 
However, from a microscopic point of view in which the plasma is modeled as 
a system of discrete charges, such an assumption has to be justified. In such 
a perspective, diamagnetism is a dynamical property that can be present 
or absent, according to the motions being ordered or chaotic. In fact, exis- 
tence of magnetization corresponds to the prevailing of gyrational motions 
transverse to the field, whereas in the state of statistical equilibrium (i.e.. 
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with prevailing chaotic motions) magnetization vanishes (see for example 
|21j). This breakdown is thus a global characteristic feature of magnetized 
plasmas, irrespective of the particular mechanism employed for obtaining 
confinement. 

The conception that magnetization due to orbital motions can exist only 
in a nonequilibrium state, characterized by motions of ordered type, was 
apparently first proposed by Bohr (see [22j, page 382). Now, Bohr took for 
granted that the relaxation time to equilibrium would be very short, as "the 
collective motions of the electrons would disappear very rapidly". On the 
other hand, we are well acquainted with the fact that the relaxation time 
from order to chaos can be very long, as occurs for example with glasses and 
with the FPU model, and was recently pointed out also in connection with 
orbital magnetization [23]. Thus, in order to establish up to which time is the 
magnetized state conserved, one should estimate a typical relaxation time 
after which the system becomes chaotic (see for example the "characteristic 
time of mixing" defined in [21], sect. 5). 

It is well known that this is a quite hard task. However, estimates of 
the relaxation time from below are available through perturbation theory, 
as we now recall. Indeed, in general such a theory allows one to construct 
adiabatic invariants /^"^ at any order n, providing for their changes I^"'^ — /(") 
(where Xt denotes the time evolved at time t of any dynamical variable X) 
estimates which in their simplest form are of the type 



_ j(n) 



<n!e"+iZ-, (1) 



where e is the perturbation parameter, while r and I are a characteristic time 



j{n) _ j{n) 



< 



and a characteristic value of I, of the system. Now, imposing 

le, and recalling n! ~ {n/eY^, formula ([1]) gives t < r(e/ne)" for all n, which, 
by taking the optimal value of n, n(e) ~ 1/e, gives t < rexp(l/e). Thus, a 
lower estimate to the relaxation time is obtained, which is exponentially long 
in 1/e as long as e < 1. It is thus clear that the condition e = 1 provides a 
natural stochasticity threshold, which should identify the relevant transition, 
at least as concerns the order of magnitude of the characteristic parameters 
of the problem. Indeed, for smaller e the motions keep an ordered character 
for practically infinite times, whereas for larger e the ordered character is 
not even guaranteed up to the microscopic time r. As a matter of fact, the 
estimates for the changes of the adiabatic invariants are in general a little 
more complicated than ([T|), and the lower estimates for the relaxation time 
are found to increase as stretched (rather than pure) exponentials, but the 
conclusion for the stochasticity threshold to be drawn in a moment, remains 
unaltered. 

This classical scheme was implemented in a probabilistic frame in the 
paper [T6]. Later the scheme was shown to be applicable also for systems 
of macroscopic sizes, i.e., in the so called thermodynamic limit [T7 1 [T8l [TU] , 
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which is an essential point for our purposes. In such a probabihstic frame, 
one renounces to control the changes of the adiabatic invariant along all 
single trajectories, and just controls mean properties with respect to a given 
invariant measure in phase space. For example, one can look at the time au- 
tocorrelation function Cj{n) (t) of the adiabatic invariant at order n, defined 
as usual by 



where (•) denotes mean with respect to the given measure. In terms of the 
time autocorrelation function, the analogue of the classical estimate ([1]) then 
takes the form 



where = (-^^) — i^)'^ is the variance of X. The latter provides a natural 
dimensional constant for the autocorrelation, since one has Cx(0) = CJ^. By 
optimization with respect to n, the time after which the adiabatic invariant 
may lose correlation is still found to be exponentially long in 1/e, provided 
one has e < 1. So, e = 1 again turns out to be the perturbation estimate of 
the stochasticity threshold. 

Our main task is thus to estimate the stochasticity threshold for a mag- 
netized plasma, in the probabilistic frame just sketched. In section [3] we 
describe the model we study, define the dynamical variable of interest (the 
component of the angular momentum of each electron along the field) and 
give the lowest order estimate for its time autocorrelation function. This 
leads to a natural conjecture for identifying the perturbation parameter e, 
which then gives the stochasticity threshold by the condition e = 1. In sec- 
tion m we give a general formula for the threshold in terms of temperature 
and of the fluctuations of the microfield. Using the available analytical esti- 
mate of such fluctuations for the model of a one component plasma with a 
neutralizing background |20] . a deflnite formula for the theoretical density 
limit is then obtained. In section [S] the theoretical density limit is compared 
to the empirical data for collapses in fusion machines. Some comments are 
flnally added in the conclusions. 

3 Definition of the model. Conjecture on the stochas- 
ticity threshold 

The model chosen for the magnetized plasma is the simplest one we could 
conceive in order to check the main idea of the present paper, namely, that 
the relevant feature concerning loss of conflnement, regardless of the par- 
ticular mechanism involved in each machine, is the occurring of a sharp 
transition from order to chaos as the perturbation due to the microfield is 



c,(„)(t) = (4")/("))-(/W)2 
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increased beyond a threshold. So, first of all, for what concerns the mag- 
netic field B, we consider the extremely idealized case in which it is uniform, 
say B = Bez, where is the unit vector along the z axis. Concerning the 
plasma itself, the key point is that it should be conceived as a dynamical 
system of point charges, and not as a continuum. Thus, any charge will be 
subiect, in addition to the Lorentz confining force due to B, also to the force 
of the microfield E, defined as the vector sum of the Coulomb fields created 
by all the other charges. Mutual magnetic forces and retardation effects are 
neglected. 

So we have a system of several kinds of charges, and the Newton equation 
for the jth charge (in the nonrelativistic approximation) is then 

^j^j = ejvj A B + ejEj (2) 

where mj and Cj are the mass and the charge of the particle, Xj and Vj = icj 
its position vector and velocity, and Ej = E(xj) the microfield (evaluated 
at Xj, which depends on the positions of all the charges. Finally, in order 
that the dynamical system be defined within the standard approach of er- 
godic theory, we consider as given also an invariant measure, a few minimal 
properties of which will be mentioned later. 

If the microfield is neglected, the transverse motion of each particle is a 
uniform gyration about a field line with its characteristic cyclotron frequebcy 
Ucj = \ej \ B/ruj. So the system is integrable, the z component of the angular 
momentum of each particle being a constant of motion. The microfield, 
acting as a perturbation, makes the system no more integrable. 

For what concerns the adiabatic invariant to be investigated, in principle 
we should look at the magnetization of the system, to which each charge 
contributes through the z component of its angular momentum. However, 
it is well known that only the electrons are relevant, the ions contribution 
to magnetization being negligible. So we will consider the contribution to 
magnetization due to any single electron, i.e., the z component of its angular 
momentum. Since now on, the index j referring to a chosen electron will 
be left understood. Thus, as zeroth order approximation for the adiabatic 
invariant we take the quantity 



(v^ denoting transverse velocity of the chosen electron), which is propor- 
tional to the transverse kinetic energy of the electron. One immediately 
checks (see page 16 of the book of Alfven [23], or any plasma physics text- 
book) that L is the z component of the angular momentum of the chosen 
electron, referred to its instantaneous gyration center (or guiding center), 
the latter being calculated in the approximation in which the perturbing 
force is neglected 
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We also add here the formula for the time derivative L of L, as we will 
need it in a moment. As L is a multiple of • v^, L is immediatley obtained 
through dot multiplication of Newton's equation ([2]) by 2v_l, which gives 

• 2m „ 

L = — vx • Ex . (4) 

We come now to the main point: to find the dimensionless perturbation 
parameter e which determines the stochasticity threshold corresponding to 
the destruction of the chosen adiabatic invariant (and of all the adiabatic 
invariants corresponding to each electron). This would require performing 
the corresponding perturbation estimates at all orders, which at the moment 
we are unable to do. What we can easily do is to perform the zeroth order 
estimate for the time autocorrelation function of L, which turns out to be 



cut) ^ , 1 



a 
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4 s 1-5^ P) 

Notice that as characteristic microscopic time r of the unperturbed elec- 
tron's motion we have naturally taken l/wc- 

The proof of ([5]) is rather simple. One starts from the elementary identity 

CL{t)=al-\{{Lt-Lf) , 

and uses the inequality {{Lt — Vf') < (L^)t^, which is just a function the- 
oretic analogue of the Lagrange finite increment formula of elementary cal- 
culus, and basically follows from unitarity of the time evolution of the dy- 
namical variables (see [T7], Theorem 1, and p!B], Section 7). This already 
gives inequality ([5]), with (L^) in place of cr^ . Relation (l5|) then immediatley 
follows by noting that, due to the time-invariance of the measure, for any 
dynamical variable X one has {X) = 0, so that (X^) = cj^. 

As explained in the Introduction, from ([5]) we are led to conjecture that 
the relevant dimensionless parameter of the problem is 

« = ^, (6) 

{o'x = denoting standard deviation), and this leads to a stochasticity 

threshold given by 

^ = 1 . (7) 

We add now a comment of a general character, which concerns the way 
in which a relaxation time for L (or analogously for any variable X) turns 
out to be identified in the present approach, which combines perturbation 
and statistical mechanics methods. Prom ([5]) one sees that the relaxation 
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time tl^^ of L is given by t^*^' = a^/a^, a formula which involves standard 
deviations. Now, compare such a formula with the one generally met in text- 
books, i.e., = L/L, where it should be understood that "typical values" 
are to be taken for the numerator and the denominator. But this requires a 
great ingenuity from the part of the reader, especially when variables are in- 
volved which have vanishing mean. So one may say that the identification of 
the relaxation time provided by perturbation theory in a probabilistic frame, 
namely, t£^' = ai/cr^, appears to be some definite quantitative implemen- 
tation of the intuitive idea underlying the familiar informal definition, and 
amounts to the prescription that the informal qualification "typical values" 
should be understood in the sense of "standard deviations" . 

We add now a final remark in which the previous comment is used to read 
in a quite transparent way the condition ([7|) which defines the stochasticity 
threshold. Indeed, through formula Q of the next section it will be seen that 
the condition for the threshold can be expressed in the form {ge^/ Ba^) = 1. 
Thus, just in virtue of the previous comment relating standard deviations 
and typical values, one sees that the threshold occurs when the typical value 
of the perturbing force due to the microfield equals the typical value of the 
Lorentz force which characterizes the unperturbed motions. So the condition 
e = 1, which we have assumend as a definition of the threshold within a 
rather abstract point of vue, is just what one would immediately guess, as 
the naivest implementation of the idea that a threshold occurs when the 
perturbing force equals the unperturbed one. 

4 The stochasticity threshold in terms of macro- 
scopic parameters. The theoretical density limit 

Our aim is now to express the stochasticity threshold ([7]) in terms of the 
macroscopic parameters T, n, temperature, electron number density and 
field strength. Recalling the expressions (l3|), ^ of L and L, and the defini- 
tion (jjc = |e| B/m, the threshold ([7|) takes the form 

B cr^2^ 

It is clear that the standard deviations appearing in ([5]) depend on the 
model of plasma adopted, which determines the microfield, as well as on the 
chosen invariant measure. The choice of the invariant measure is a quite 
delicate problem, particularly in a nonequilibriun situation as the one we 
are discussing here. A general introduction may be found in the book |26j . 
For example, it is obvious that cJvx ^'^^ should be expressed in terms 
of temperature, albeit with coefficients which depend on the assumptions 
made for the velocity distribution. Analogously, the statistical properties of 
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Figure 1: Density limit values vs B for various machines: conventional 
tokamaks, for which recent data are shown (see references [291 EOl EH [32l 
[MllSllMllMlETllMllMlllQlillllSlllallll]) along with the original ones 
of Murakami (see [2E]), stellar ator machines [l5l |36l HH HH] , and spherical 
tokamaks [50l [51] . Dotted line is the theoretical density limit (jlip . 



the microfield may be different for a system composed by electrons plus a 
neutralizing background, rather than for a system of electrons and ions. 

Quite natural assumptions on the measure are: i) that velocities and 
positions are independent variables; ii) that the distribution of the transverse 
velocities is Maxwellian at a temperature T; iii) that the distribution of 
positions is isotropic. Under these natural assumptions the equation ([8]) for 
the threshold is seen to take the form 



3 ByJkBT/m 
Indeed, from i) and ii) one gets 



1 (9) 



^2 ^2 



the variance of a vector F being defined by Cp = Gp^ + Up^ + ap^ . One also 



gets 

1 2 _ ^bT 

2 ""^^ " m 
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{ks being the Boltzmann constant) and furthermore, as one easily checks, 



Finally, from iii) one gets = (2/3) cr^ 

The form Q of the equation for the thresold already constitutes in our 
opinion a significant result. Indeed, the fluctuation of the microfield 
should in principle be itself a measurable quantity, which depends on the 
macroscopic state of the plasma, namely, electron number density n and the 
temperatures of the several constituents. So the previous relation provides 
in principle the density limit as a function of the macroscopic state of the 
plasma. 

However, we were unable to find in the literature sufficient experimental 
information on the fluctuation of the microfield. So, in order to have 
a definite theoretical formula to be compared with the experimental data, 
we limit ourselves to the consideration of a particular model for which an 
estimate of is available. In fact a formula for was given by Iglesias, 
Lebowitz and MacGowan [20j for the model of a one component plasma with 
neutralizing background, at equilibrium with respect to the Gibbs distribu- 
tion, for which they found 



where eo is the vacuum dielectric constant and n the electron number density 
(see [20], formula (2.5), substituting n for p and I/sq for An). 

So, for a one component plasma with neutralizing background at Gibbs 
equilibrium at temperature T, the stochasticity threshold (0) takes the form 



in which temperature disappeared, so that the threshold only involves den- 
sity and field strength. Notice however that this might not be true with 
a more realistic model of a plasma, in which the temperature appearing 
in pop, which refers to the plasma as a whole, may be different from the 
electron transverse temperature which enters the previous formulas. 

Formula (jlip for the limit density (holding for a one component plasma 
with neutralizing background, at Gibbs equilibrium at temperature T) is 
the type of result we were looking for, inasmuch as it provides a definite 
theoretical formula for the density limit that can be compared to the avail- 
able empirical data for collapses in fusion machines, as will be done in the 
next section. 
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5 Comparison with the empirical data for plasma 
collapses in fusion machines 



We now check whether the transition from order to chaos discussed here has 
anything to do with the empirical data for collapses in fusion machines. We 
recall that a proportionality of the density limit to the square of the magnetic 
field in tokamaks was suggested by Granetz [27] on the basis of empirical 
data, but apparently was not confirmed by later observations \4:4\ llOj . It 
is well known that, while at first a proportionality to the magnetic field 
(through B / R, where R is the major radius of the torus) had been proposed 
on an empirical basis for tokamaks by Murakami [28] , in the plasma physics 
community the common opinion is rather that the density limit for tokamaks 
should be proportional to the Greenwald parameter /p/r^, where Ip is the 
plasma current and ra the minor radius of the torus (see [TU] ) . 

We do not enter here a discussion of this point, and only content our- 
selves with plotting in Figure [T] a collection of available data of the density 
limit for several fusion machines versus their operating magnetic field B 
in log-log scale, comparing the data to the theoretical formula (jlip . The 
first thing that comes out from the figure is that the order of magnitude 
of the theoretical threshold is correct, and this without having introduced 
any phenomenological parameter. There is no adjustable parameter in the 
theory, and no fittind at all. One is thus tempted to say that the essence of 
the phenomenon has perhaps been captured, especially in consideration of 
the extreme simplicity of the model, with respect to the variety of machines 
and of operational conditions to which the experimental data refer. 

Entering now in some more details, one sees that the theoretical law 
appears to correspond not so badly to the data for the high field machines 
(tokamak and stellarators), whereas a sensible discrepancy is met for the 
low field machines (spherical tokamaks), for which the experimental data 
are larger by even an order of magnitude. Perhaps this discrepancy might 
be attributed to the fact that we are discussing here a model describing an 
isolated, non sustained, system (i.e., with no input heating power), whereas 
the low field machines considered in the figure are just the ones character- 
ized, in general, by lower confinement time and thus by larger sustainment. 
Indeed, (see the empirical Sudo limit for stellarators [38]) larger densities 
are expected to be accessible as the input power is increased (although this 
is not so clear for tokamaks [lO]). This is illustrated, in the figure, by the 
three points reported for the same device (the stellarator WS-A7 [IS]) at es- 
sentially the same applied field, which however correspond to three different 
(increasing) input heatings. 



11 



6 Conclusions 



In view of the lack of any first principles rationale for the existence of a 
density limit in fusion machines, the comparison between theory and exper- 
iments exhibited in Figure [1] appears encouraging. Particularly so, if one 
considers the extreme simplicity of the model (uniform plasma in a uniform 
field), with respect to the the variety of machines and of operational condi- 
tions to which the experimental data refer. The essence of the phenomenon 
seems to have been captured. 

So, one might consider as plausible the main proposal advanced in the 
present paper, namely, that the density limit characterizing the empirical 
collapses of fusion machines corresponds to a transition from order to chaos 
in the following sense. At low densities, ordered motions due to the imposed 
magnetic field prevail, with the electrons performing transverse gyrational 
motions, and thus with a magnetic pressure. Then, as density is increased, 
the perturbations caused by the fluctuations of the microfield (which increase 
as the density) introduce a chaotization, until a stochasticity limit (and so a 
density limit) is attained, beyond which ordered motions are lost, together 
with magnetic pressure and confinement. 

A key feature of the present approach, with respect to treatments involv- 
ing the continuum approximation, such as magnetohydrodynamics, is that 
we are dealing here with the plasma as a discrete system of charges. Indeed 
in our treatment an essential role is played by the microfield acting on a 
single electron, and so it is not clear how the instability found here could 
find place within the continuum approximation, or any other approximation 
involving high-frequency cutoffs. For an analogous role of discreteness of 
matter in cosmology, see [52] and [53] . 

Actually, even in plasma physics theory there exists a huge literature 
in which the discrete nature of matter is taken into account, following the 
approach of kinetic theory (see for example [54]). A comparison with the 
results obtained here within the approach of dynamical systems theory would 
thus be in order. We hope to come back to this problem in the future. 

A further remark is that the existence of a density limit proportional to 
the square of the magnetic field is well known in the frame of nonneutral 
plasmas (see [55]), under the name of Brillouin limit. The physical context 
is however rather different, because the density limit in the latter case refers 
to the existence of a particular motion, in which the plasma, dealt with as a 
continuum, performs a rigid rotation about the z axis. Actually, it is clear 
that a magnetization threshold in the sense discussed here should exist for 
nonneutral plasmas too. The only problem is that we are unaware of any 
estimate of the standard deviation of the microfield in such a case. We hope 
to come back to this problem in the future. 

We finally point out that the proportionality of the density limit to 
the square of the magnetic field predicted by the theoretical law ([IT]) , if 
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confirmed, might have relevant implications for future tokamaks. 
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